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(57) Abstract: The invention relates to a method for the interactive simulation 
of contact between objects. The inventive method comprises the following steps, 
namely: the parameters describing the physical characteristics of each of the objects 
are calculated; at the beginning of each simulated model sampling time period, each 
object is subjected to a real-time analysis of the specific behaviour thereof accord- 
ing to a free movement that does not take account of possible subsequent contacts, 
and, subsequently, at a global scene level, pairs of detected intersecting objects are 
subjected to a real-time analysis; a list of collision groups is established; for each 
collision group, parameters representing the physical characteristics of the objects 
and the description of the collisions are returned in real time, such as to characterise 
the contact between two objects in the case of a pure relative sliding movement; and, 
for each object, the specific behaviour of the object following the collision is dis- 
played in real time and the set of real-time processes is performed with a shorter 
calculation time than the sampling time. 

(57) Abrege : Selon le precede de simulation interactive du contact entre objets on 
calcule les parametres decrivant les caracteristiques physiques de chacun des ob- 
jets ; au debut de chaque pas de temps d'echantillonnage d'un modele simule, on 
precede au niveau de chaque objet a une analyse en temps reel du comportement 
propre de l'objet selon un mouvement libre qui ne tient pas compte d'eventuels 
contacts ulterieurs, puis on analyse en temps reel, au niveau d'une scene globale, 
des paires d'objets qui sont detectes en intersection ; on etablit une liste de groupes 
de collisions ; on rapatrie en temps reel, pour chaque groupe de collisions, des pa- 
rametres representant les caracteristiques physiques des objets et la description des 
collisions, de maniere a caracteriser le contact entre deux objets dans le cas d'un 
glissement relatif pur ; on precede au niveau de chaque objet a une visualisation en 
temps reel du comportement propre de l'objet suite a la collision, et 1 'ensemble des 
traitements en temps reel s'effectue avec un pas de temps de calcul plus court que 
le pas de temps d'echantillonnage. 
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Precede et dispositif de simulation interactive du contact entre obiets 

La presente invention concerne un procede et un dispositif de 
simulation interactive du contact entre au moins un premier objet 
5 deformable et au moins un second objet avec un pas predetermine de 
temps d'echantillonnage d'un modele simule. 

On a deja propose d'effectuer une simulation de mesures 
d'interpenetration entre un objet rigide et un objet deformable a partir 
d'estimations volumiques ou de distance, notamment pour des 
10 applications de chirurgie virtuelle ou un outil virtuel rigide de chirurgie 
coopere avec un organe virtuel deformable du corps humain. 

Toutefois, selon ces methodes la relation entre la mesure 
d'interpenetration et les forces de contact de reaction n'ont pas de base 
physique et des forces artificielles peuvent s'appliquer a des nceuds des 
15 maillages des objets qui ne sont pas en contact, ce qui nuit a la fiabilite, 
ou les forces de contact ne respondent pas aux conditions du probleme de 
Signorini. 

La presente invention vise a remedier aux inconvenients precites et 
a permettre de realiser une simulation interactive en temps reel du contact 

20 entre des objets, dont au moins certains sont deformables, de facon 
simplifiee et economique tout en respectant les contraintes des lois 
physiques qui regissent les contacts, de telle maniere que les contacts 
simules entre des objets soient fiables et qu'ainsi la stabilite de la 
simulation soit garantie. 

25 Ces buts sont atteints grace a un procede de simulation interactive 

du contact entre au moins un premier objet deformable et au moins un 
second objet avec un pas predetermine de temps d'echantillonnage d'un 
modele simule, 

caracterise en ce que : 

30 (a) on calcule au prealable les parametres decrivant les caracteristiques 
physiques de chacun des objets, tels que la geometrie et la 
mecanique des materiaux de chacun des objets, et on stocke ces 
parametres dans une memoire, 

(b) au debut de chaque pas de temps d'echantillonnage du modele 
35 simule, on procede au niveau de chaque objet a une analyse en 
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temps reel du comportement propre de I'objet pour predire les 
positions, vitesses et accelerations de cet objet selon un mouvement 
libre qui ne tient pas compte d'eventuels contacts ulterieurs, 

(c) a chaque pas de temps d'echantillonnage du modele simule, on 
5 analyse en temps reel, au niveau d'une scene globale comportant les 
objets susceptibles de venir en contact, des paires d'objets qui sont 
detectes en intersection, et on etablit une liste de groupes de 
collisions qui contient une chaTne d'objets en collision et une 
description des collisions, 

10 (d) a chaque pas de temps d'echantillonnage du modele simule, on 
rapatrie en temps reel, pour chaque groupe de collisions, des 
parametres representant les caracteristiques physiques des objets et 
la description des collisions, de maniere a determiner, pour chaque 
cas, la solution au probleme de Signorini qui gouverne le contact 

15 entre deux objets dans le cas d'un glissement relatif pur, 

(e) a la fin de chaque pas de temps d'echantillonnage du modele simule, 
on procede au niveau de chaque objet a une visualisation en temps 
reel du comportement propre de I'objet suite a la collision, et 

(f) I'ensemble des traitements en temps reel s'effectue avec un pas de 
20 temps de calcul plus court que le pas de temps d'echantillonnage du 

modele simule, de maniere a definir une simulation interactive ou 
I'utilisateur peut intervenir directement en cours de simulation. 

Lors de I'etape a) de calcul prealable des parametres decrivant les 
caracteristiques physiques de chacun des objets, on utilise pour les 
25 parametres decrivant la mecanique des materiaux une description des 
deformations de type elements finis, avec le remplissage et I'inversion de 
matrices, la resolution de systemes d'equation et le stockage des donnees 
en memoire. 

Selon un mode particulier de realisation, chaque objet est decrit 
30 dans une configuration au repos comme un ensemble de triangles 
reproduisant sa surface et un ensemble de tetraedres decrivant I'interieur 
de I'objet. 
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Avantageusement, chaque triangle est decrit par trois points, places 
dans un ordre qui permet de calculer des normales qui sont 
invariablement dirigees vers I'exterieur de I'objet. 

De preference, les deformations des objets sont interpolees par la 

5 methode des elements finis en utilisant un maillage tetraedrique lineaire. 

A chaque pas de temps de calcul, on integre lors de I'etape b) au 
niveau d'un objet, les forces explicites appliquees sur I'objet, qui sont deja 
connues au depart du pas de calcul, de maniere a definir le mouvement 
qu'elles creent sur I'objet, tandis que la valeur des forces de contact 

10 implicites, qui dependent elles-memes du mouvement des objets dans le 
pas de temps de calcul, est determinee lors de I'etape d) de recherche au 
niveau d'une scene globale, de la solution au probleme de Signorini. 

Lors de I'etape c) d'analyse au niveau d'une scene globale, on detecte 
geometriquement les intersections existantes entre les objets de la scene 

15 pour extraire des couples d'elements d'objets en intersection, une 
longueur et une direction d'interpenetration entre les deux elements d'un 
couple d'elements d'objets. 

Selon une variante de realisation, lors de I'etape c) d'analyse au 
niveau d'une scene globale, pour extraire des couples d'elements d'objets 

20 en intersection, une longueur et une direction d'interpenetration entre les 
deux elements d'un couple d'elements d'objet, on prend egalement en 
compte un mouvement intermediate des objets entre le pas de calcul 
precedent et le pas de calcul courant, pour calculer une direction 
d'interference privilegiee entre ces objets. 

25 Lors de I'etape d) de recherche de la solution au probleme de 

Signorini, on reconstruit les points extremes d'application de la force de 
contact entre deux objets soumis a une collision lorsque ces points 
extremes d'application n'ont pas ete determines a I'etape precedente. 

Selon un mode de realisation particulier, lors de I'etape d) lors de 

30 I'etape d) de recherche de la solution au probleme de Signorini, dans le 
cas d'une intersection segment-segment de deux objets en triangle, les 
deux points choisis pour constituer les points extremes d'application de la 
force de contact entre les deux objets soumis a une collision sont situes a 
I'intersection de chacun des deux segments avec le plan forme par la face 

35 du triangle en intersection. 
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Selon un autre mode de realisation particulier, lors de I'etape d) de 
recherche de la solution au probleme de Signorini, dans le cas d'une 
intersection point-face de deux objets en triangle, un premier point choisi 
pour constituer un point extreme d'application de la force de contact entre 
les deux objets soumis a une collision est le point de I'intersection tandis 
que le second point extreme d'application de la force de contact entre les 
deux objets soumis a une collision est la projection du premier point 
extreme sur la face du triangle en intersection. 

Selon un aspect particulier de la presente invention, on utilise les 
coordonnees barycentriques pour repartir les deplacements et les forces 
des points d'application de la force de contact entre les points extremes 
d'application de la force de contact en effectuant une interpolation lineaire 
pour une moderation en elements finis. 

De facon plus particuliere, on peut calculer la distance 6 
d'interpenetration entre les deux points extremes d'application de la force 
de contact dans le cas d'un contact segment-segment entre un premier 
segment et un second segment d'un second triangle a partir de I'equation 
suivante : 



ou : a et 1-a sont les coordonnees barycentriques sur le premier 
segment, 

(3 et 1-p sont les coordonnees barycentriques sur le deuxieme 
segment, 

a, bi q sont les coordonnees de la direction n* d'interpenetration, 

Wi et W 2 sont les coordonnees du premier segment, 

Vi et V 2 sont les coordonnees du second segment. 

On peut egalement calculer la distance 8 d'interpenetration entre 
les deux points extremes d'application de la force de contact dans le cas 
d'un contact point-plan entre un point d'un second triangle et un plan d'un 
premier triangle a partir de I'equation suivante : 
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ou : a, (3 et y sont les coordonnees barycentriques sur le premier triangle, 

a-, bj q sont les coordonnees de la direction n, d'interpenetration, 

5 Wi, W 2 , W 3 sont les coordonnees du premier triangle, 

Vi represente les coordonnees du point de contact constitue par un 
sommet du second triangle. 

Apres avoir determine les points duplication des forces de contact 
entre deux objets en collision, lors de I'etape d) on procede au transfert 
10 des caracteristiques mecaniques des objets dans I'espace defini des 
contacts dans lequel on traite I'ensemble d'un groupe de m contacts avec 
n objets ou m et n sont des entiers. 

De facon plus particuliere, lors de I'etape d) on considere la masse 
et I'inertie d'un objet de facon globale, en son centre de gravite et on 
15 etablit une relation instantanee entre les forces de contact f c dans la 
direction du contact, les accelerations S] dues aux contraintes dans la 
meme direction et les accelerations libres 8" libre dans la meme direction 
connues lors de I'etape c) au niveau d'une scene globale, selon I'equation 
suivante : 

20 

J c M-'j c T f c+< y librc( 3) 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane lineaire et angulaire dans I'espace des contacts, 

25 3 T C est la matrice transposee de 3 C/ 

M est une matrice bloc diagonale correspondant a la masse et 
I'inertie de n objets du groupe de contacts. 

Lors de I'etape d) pour le transport des caracteristiques mecaniques 
locales, on etablit une relation entre : 
30 • la difference de deplacement (Uk') des points du maillage 

deformable representant I'objet i a I'instant k, entre la 
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deformation libre (U'wibre) et la deformation contrainte (U'k /C ) soit 

Uk = U'k,c " U'k.libre 

• Les positions relatives des objets, libres et contraintes, dans 
5 I'espace des contacts : Snbre et 5 C . 

5 = Ei=i n N c ' Uk' + 5, ibre (4) 

ou N c ' est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements aux contacts. 
10 De meme, on etablit une relation entre les forces dans I'espace des 

contacts f c et les forces dans I'espace des forces de deformation F k . 

Fk = (N c i ) T / c (5) 

De fagon plus particuliere, lors de I'etape d) on etablit une relation 
15 lineaire instantanee de caracterisation des deformations ou deplacement 
de contact 8 C a partir des forces de contact f c et des deplacements libres 
8iibre dus aux mouvements libres integrant uniquement les forces connues 
de fagon explicite en debut de pas de temps de calcul, selon I'equation 
suivante : 

20 <? = IX; n; a (u k ,> (n;> t ] / c+ \ bre (6) 

ou : est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements de contacts, 

(N^) T est la matrice transposed de N' c/ 

25 A est une matrice permettant de definir la deformation de I'objet au 

niveau local, de telle sorte que si Uk represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et Uk-i 
represente le vecteur de deplacement dans le repere local de I'objet 
au pas de calcul precedent dont les valeurs instantanees sont 

30 connues au debut du pas de calcul courant, on a : 

U K = A (Uk-i) F k + b (U k -i) (7) 



35 



ou Fk est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, et b est un vecteur qui a une 
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valeur dans I'espace des deplacements et qui depend du modele de 
deformation de I'objet. 

Dans un cas plus general, tors de I'etape d) on etablit une relation 
instantanee de caracterisation des deformations ou deplacements de 
5 contact 5 C a partir des forces de contact f c et des deplacements libres 5iibre 
dus aux mouvements libres integrant uniquement les forces connues de 
fagon explicite en debut de pas de temps de calcul, selon I'equation 
suivante : 

10 S m = [6 dt 2 J c M 1 J c T +2 i=1 " N c ' A (lU (N C ') T ] f c+ S Uhre (8) 

ou : J c est une matrice jacobienne m*6n qui transfere !e mouvement 
instantane lineaire et angulaire dans I'espace des contacts, 

3 T C est la matrice transposee de J c ' 

15 M est une matrice bloc diagonale correspondant a la masse et 

I'inertie des n objets du groupe de contacts, 

0 est une constante dependant de la methode d'integration en 
temps, 

est une matrice de passage de I'espace des deplacements du 
20 maillage vers I'espace des deplacements de contacts, 

(N^) T est la matrice transposee de N'., 

A est une matrice permettant de definir la deformation de I'objet au 
niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et U k -i 
25 represente le vecteur de deplacement dans le repere local de I'objet 

au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 

U K = A (U k -i) F k + b (U k -i) (7) 

30 

ou F k est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, et 
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b est un vecteur qui a une valeur dans I'espace des deplacements et 
qui depend du modele de deformation de I'objet. 

On peut ainsi remarquer que dans le cas d'un objet rigide 
indeformable, on a Uk = Uk-i qui traduit I'absence de modification en 
5 fonction du temps du vecteur Uk. 

Avantageusement, le procede selon invention comprend en outre 
une etape de couplage avec un module d'interface haptique pour produire 
un retour de sensation haptique sur un dispositif mecanique par lequel un 
operateur manipule les objets dans une scene virtuelle. 

10 L'invention concerne egalement un dispositif de simulation 

interactive du contact entre au moins un premier objet deformable et au 
moins un second objet avec un pas predetermine de temps 
d'echantillonnage d'un modele simule, caracterise en ce qu'il comprend : 

(a) un module de calcul prealable des parametres decrivant les 
15 caracteristiques physiques de chacun des objets, tels que la 

geometrie et la mecanique des materiaux de chacun des objets, 

(b) une memoire de stockage des parametres prealablement calcules 
dans le module de calcul, 

(c) un module de couplage avec une interface utilisateur comprenant un 
20 dispositif mecanique tenu par un utilisateur lui permettant d'exercer 

virtuellement des forces sur lesdits objets dans une scene du modele 
simule, 

(d) un ecran de visualisation pour visualiser lesdits objets representes 
sous forme de maillages, 

25 (e) une unite centrale de traitement associee a des moyens d'entree, 
comprenant au moins 

el) un module d'analyse d'objet pour analyser en temps reel au 
niveau de chaque objet le comportement propre de I'objet pour 
predire les positions, vitesses et accelerations de cet objet selon 
30 un mouvement libre qui ne tient pas compte d'eventuels 

contacts ulterieurs, 

e2) un modele d'analyse d'une scene globale comportant les objets 
susceptibles de venir en contact, pour analyser en temps reel 
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des paires d'objets qui sont detectes en interaction et etablir 
une liste de groupes de collisions qui contient une chaine 
d'objets en collision et une description des collisions, 

e3) un module de rapatriement en temps reel, pour chaque groupe 
5 de collisions, des parametres representant les caracteristiques 

physiques des objets et la description des collisions pour 
determiner, pour chaque cas, la solution au probleme de 
Signorini qui gouverne le contact entre deux objets dans le cas 
d'un glissement relatif pur, 

10 e4) un module de traitement de chaque objet pour permettre en 

temps reel au niveau de chaque objet une visualisation en 
temps reel du comportement propre de I'objet suite a une 
collision, et 

e5) des moyens de determination d'un pas de calcul plus court que 
15 le pas de temps d'echantillonnage du modele simule de maniere 

a definir une simulation interactive. 

Avantageusement, le dispositif comprend des moyens de production 
d'un retour de sensation haptique sur I'interface utilisateur. 

Selon une caracteristique avantageuse, le pas de calcul correspond 
20 a une frequence egale ou superieure a environ 500 Hz. 

D'autres caracteristiques et avantages ressortiront de la description 
suivante de modes particuliers de I'invention, donnes a titre d'exemples, 
en reference aux dessins annexes, sur lesquels : 

- la Figure 1 est un schema montrant les differentes etapes d'un 
25 procede de simulation interactive du contact entre objets, selon 

I'invention, 

- la Figure 2 est un schema montrant differents niveaux de 
traitement de I'interaction entre objets au cours de differentes etapes du 
procede de simulation de la Figure 1, 

30 - les Figures 3A a 3C represented trois exemples d'interaction entre 

deux objets virtuels representes par des triangles, 

- la Figure 4 represente une interaction entre deux objets en 
triangle dans le cas d'une intersection segment-segment, 

- la Figure 5 represente une interaction entre deux objets en 
35 triangle dans le cas d'une intersection point-face, 



WO 2005/093610 



PCT/FR2005/000699 



- 10- 

- la Figure 6 illustre schematiquement le cas d'une collision entre 
deux objets pour laquelle la configuration peut §tre definie a partir de 
seuls criteres geometriques, 

- la Figure 7 illustre schematiquement le cas d'une collision entre 
5 deux objets pour laquelle la configuration est definie en tenant compte 

d'un mouvement intermediate, 

- la Figure 8 est un schema-bloc montrant les constituants de base 
d'un dispositif de simulation interactive du contact entre objets, selon 
I'invention, 

10 - la Figure 9 montre un exemple de contact entre un objet virtuel 

deformable et un autre objet virtuel, et 

- les Figures 10A a 10C montrent trois positions relatives differentes 
entre un objet virtuel deformable, en forme de pince, et un objet virtuel 
rigide lors d'un processus de mise en place de Pobjet virtuel deformable en 

15 forme de pince sur I'objet virtuel rigide. 

La Figure 8 illustre de fagon schematique un exemple de dispositif 

permettant de mettre en ceuvre I'invention et de realiser la simulation 

interactive en temps reel du contact entre objets tout en permettant 

notamment d'avoir un retour de sensation haptique. 
20 Une unite centrale de traitement 100, qui peut etre constitute a 

partir d'un ordinateur de type classique, permet d'effectuer les differents 

calculs necessaires pour realiser une simulation. 

Un ecran de visualisation 107 relie a I'ordinateur 100 par une 

interface graphique permet 1'affichage d'objets representes sous la forme 
25 de maillage comprenant des nceuds ou sommets reliant des segments ou 

aretes. 

Des informations peuvent etre fournies a I'ordinateur 100 a partir 
d'une interface utilisateur classique 103 pouvant comprendre un clavier et 
par exemple une souris et constituant des moyens d'entree. 

30 Un dispositif mecanique specifique 104 tenu par un utilisateur relie 

par un module de couplage 101 a I'ordinateur 100 peut en outre etre 
prevu pour permettre a I'utilisateur d'exercer virtuellement des forces sur 
les objets dans une scene d'un modele simule. Un tel dispositif mecanique 
104 et le module de couplage 101 constituent une interface haptique qui 

35 permet a I'utilisateur d'exercer une sollicitation sur les objets virtuels de la 
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scene et de recevoir en retour une simulation haptique qui est une 
reponse fournie par la simulation du contact entre objets. 

Uordinateur 100 comprend de fagon classiquement au moins un 
processeur, une memoire permanente de stockage de programmes et de 
5 donnees et une memoire de travail cooperant avec le processeur. Des 
supports externes de memoire (disquettes, CD-ROM,...) ou un modem de 
liaison a un reseau peuvent naturellement etre utilises pour charger dans 
Pordinateur des programmes ou donnees permettant d'effectuer tout ou 
partie du traitement de simulation. Sur la Figure 8, on a simplement 

10 represents a titre symbolique un exemple de memoire de stockage 102 
cooperant avec le module 100 et pouvant etre constitute par Tun ou 
Pautre des types de memoires sus-mentionnes. 

D'une maniere generate, au debut d'un processus de simulation les 
parametres decrivant la geometrie et la mecanique des materiaux des 

15 objets a simuler sont calcules dans Punite centrale 100 et stockes dans 
une zone memoire de la memoire 102. 

Pour caracteriser les deformations mecaniques des objets, il est 
utilise lors du traitement par Punite centrale 100 une description des 
deformations de type elements finis. Cela se traduit par le remplissage et 

20 Tinversion de matrices, la resolution de systemes d'equations et le 
stockage des donnees dans la memoire 102 associee a Punite centrale 
100. 

Les positions et formes courantes des objets sont evaluees en 
fonction des solicitations exercees et des lois mecaniques qui regissent les 

25 objets de la scene. 

Selon Pinvention, pour garantir une simulation stable, il est pris en 
compte dans le calcul des contacts simules entre objets les lois physiques 
qui regissent le contact. Pour permettre une simulation en temps reel, 
c'est-a-dire avec un delai tres court et borne separant la sollicitation 

30 exercee par Putilisateur via ('interface haptique 104, 101 et la reponse 
fournie par Punite centrale de traitement 100 a cette interface haptique 
104, 101, le dispositif de simulation met en oeuvre trois modules 
principaux qui sont sollicites iterativement a chaque pas de temps 
d'echantillonnage du modele simule. Par ailleurs, Pensemble des 

35 traitements en temps reel s'effectue avec un pas de temps de calcul plus 
court que le pas de temps d'echantillonnage du modele simule. 
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Les trois modules principaux du dispositif de simulation mettant en 
oeuvre les differentes etapes du procede de simulation se presentent 
essentiellement de la facon suivante : 

un premier module dit "mecanique" se situant au niveau de chaque 
5 objet et decrivant son comportement propre, permet de faire evoluer 

la position et la forme de I'objet suivant les forces et les lieux des 
forces exerces. Ce module est appele en debut de pas de calcul pour 
predire les positions, vitesses et accelerations des objets sans tenir 
compte du contact puis sera de nouveau mobilise pour tenir compte 
10 des forces calculees dans un troisieme module dit "de traitement du 

contact". 

Un second module dit de "detection de collision", se situant au 
niveau de la scene globale etablit des paires d'objets qui sont 
detectes en intersection. Ce module, de facon optionnelle, peut creer 

15 des mouvements intermediates entre les pas de la simulation pour 

savoir quand et comment les objets sont entres en intersection. Ce 
module est regi avant tout par des lois geometriques optimisees qui 
permettent d'accelerer le calcul afin d'obtenir une chaine d'objets en 
collision et une description des collisions. Un groupe de collision est 

20 ainsi un ensemble d'objets relies entre eux par au moins une 

collision. Un objet entre dans un groupe s'il est en collision avec au 
moins un des objets du groupe. Une collision est decrite 
obligatoirement par la paire d'objets en collision et par le lieu de la 
collision en utilisant soit les elements geometriques de base (par 

25 exemple deux triangles ou deux surfaces) en intersection, soit par un 

segment reliant les deux points qui sont localement les plus 
interpenetres. 

Un troisieme module dit "traitement du contact" est appele par le 
module de "detection de collision" et fait appel en retour au module 

30 "mecanique". Pour chaque groupe de collisions, le module de 

traitement du contact rapatrie les caracteristiques physiques des 
objets et la description des collisions. Le module est apte a 
determiner, pour chaque cas, la solution au probleme de Signorini 
qui gouverne le contact entre deux objets dans le cas du glissement 

35 pur. 
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Uinvention permet d'effectuer une simulation interactive. Une 
simulation est definie par le pas de temps d'echantillonnage du modele 
simule et par le pas de temps de calcul de ce modele. Le precede selon 
Pinvention met en oeuvre un pas de temps de calcul qui est toujours 

5 inferieur au pas de temps choisi pour Pechantillonnage, ce qui permet 
d'avoir une simulation interactive ou Putilisateur va pouvoir intervenir 
directement en cours de simulation. 

La Figure 1 resume les principales etapes du precede selon 
Pinvention qui met en oeuvre une boucle de simulation utilisant les trois 

10 modules principaux precites installes dans Pordinateur 100 de la Figure 8. 
La Figure 2 illustre les differents niveaux de traitement entre objets au 
cours des differentes etapes du procede de simulation. 

Une premiere etape de traitement 130 utilise le module dit 
"mecanique" et se situe au niveau de chaque objet (niveau objet 3). Les 

15 informations sont fournies a travers un module de couplage 120 a partir 
de Pinterface utilisateur 110 ou interface haptique qui determine la 
position et la forme de chaque objet (informations 135 elaborees a Petape 
130). 

Dans cette premiere etape 130, une moderation 13 prend en 

20 compte chaque objet ou outil 201, 202, 203 de fagon individuelle sans 
tenir compte d'interactions ulterieures eventuelles et permet de faire 
evoluer la position et la forme de Pobjet suivant les forces et les lieux des 
forces exerces a partir de Pinterface utilisateur 110. 

Une deuxieme etape du traitement 140 utilise le module dit 

25 "detection de collision" et se situe au niveau d'une scene globale (niveau 
scene 4). Les informations 135 elaborees a Petape 130 sont utilisees a 
Petape 140 pour etablir des paires d'objets qui sont detectes en 
intersection. Au cours de cette etape 130, il est elabore une liste de 
groupes de collisions (informations 145) contenant une chame d'objets en 

30 collision et une description des collisions. Dans cette etape 140, une 
moderation 14 prend ainsi en compte une paire d'objets en intersection 
tels que les objets 201, 202 au niveau d'une scene globale. 

Une troisieme etape de traitement 150 utilise le module dit 
"traitement du contact" et se situe au niveau d'une scene globale (niveau 

35 scene 5). Les informations 145 elaborees a Petape 140 ainsi que les 
informations 135 elaborees a Petape 130 sont utilisees pour determiner 
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pour chaque cas la solution au probleme de Signorini qui gouverne le 
contact entre deux objets dans le cas de glissement pur (information 155). 
Dans cette etape 150, une modelisation 15 prend ainsi en compte 
I'interaction entre deux objets tels que les objets 201, 202 au niveau d'une 

5 scene globale, les caracteristiques physiques des objets et la description 
des collisions etant rapatriees pour chaque groupe de collisions par le 
module "traitement de contact". 

La troisieme etape de traitement 150 fournit des informations 155 
concernant des forces et des lieux qui sont transmises au premier module 

10 "mecanique" au cours d'une quatrieme etape de traitement 160 qui se 
situe a nouveau au niveau des objets (niveau objet 6). Dans cette etape 
160, le resultat du traitement de simulation en temps reel peut etre 
simplement normalise dans une etape 170 de visualisation ou peut §tre 
transmis en retour a travers le couplage 120 vers I'interface utilisateur 110 

15 pour conferer a I'utilisateur un retour de sensation haptique. Dans cette 
etape finale, une modelisation 16 prend ainsi a nouveau en compte 
chaque objet ou outil 201, 202, 203 individuellement tout en ayant tenu 
compte des contacts precedemment simules. 

Au niveau de I'objet, dans une forme de realisation preferee, on 

20 decrit I'objet comme un ensemble de triangles reproduisant sa surface et 
un ensemble de tetraedres pour decrire son interieur, le tout dans une 
configuration de repos. Cette configuration correspond a la forme de 
I'objet lorsqu'aucune force ne lui est appliquee. 

Avantageusement, les triangles sont decrits par trois points, places 

25 dans un ordre qui permet au calcul des normales d'etre invariablement 
dirigees vers I'exterieur de I'objet. Les surfaces des objets sont fermees de 
facon a pouvoir distinguer un exterieur d'un interieur. 

Les deformations des objets sont interpolees par la methode des 
elements finis en utilisant un maillage tetraedrique lineaire. Le dispositif 

30 permet de simuler differentes lois de comportement sous reserve que Ton 
puisse en extraire localement, approximativement et pour un pas de 
calcul, une relation lineaire entre les forces exercees et les deplacements 
autour d'une configuration locale. 

Si Ur represente le vecteur de deplacement dans le repere local 

35 d'un objet a I'instant courant t et si Uk-i represente le vecteur de 
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emplacement dans le repere local de I'objet au pas de calcul precedent t-1, 
on a la relation suivante : 

U k = A(U k -i) F k + b(U k .i) (7) 

5 

ou : A est une matrice permettant de definir la deformation de I'objet au 
niveau local, 

Fk est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, 
10 b est un vecteur qui a une valeur dans I'espace des deplacements et 

qui depend du modele de deformation de I'objet, et 
Uk-i est un vecteur dont les valeurs instantanees sont connues au 
debut du pas de calcul, de Pinstant courant t 

15 Avantageusement, une distinction est faite entre le mouvement 

global de Pobjet decrit par la relation fondamentale de la dynamique et la 
deformation de I'objet autour d'une configuration courante decrite par la 
loi de deformation. 

Les forces appliquees sur I'objet se distinguent par leur caractere 

20 explicite ou implicite. Une force explicite est d'ores et deja connue au 
depart du pas de calcul, et il faut I'integrer pour connaitre le mouvement 
qu'elle cree sur I'objet. Les forces de contact, au contraire, sont implicites 
au sens ou elles dependent elles-memes du mouvement des objets dans 
le pas de temps. Pour un pas de temps, on integre done les forces 

25 explicites et on passe alors au niveau de la scene globale pour pouvoir 
trouver la valeur des forces implicites. 

Apres Pintegration des forces explicites au niveau de chaque objet, 
on appelle position et forme "libre" la configuration des objets dans la 
scene. Cette configuration est obtenue sans Intervention des forces de 

30 contact. Ainsi, un "mouvement libre" est un mouvement integrant 
uniquement les forces connues de facjon explicite en debut de pas de 
temps de calcul. Par suite on considere ici que le mouvement libre est le 
mouvement obtenu sur un pas de temps quand les forces de contact ne 
sont pas integrees. 

35 Le systeme propose met alors en oeuvre un processus de detection 

de collision qui permet de tester geometriquement les intersections 
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existantes entre les objets 223, 230 de la scene et les directions 
privilegiees pour sortir les objets de cette collision (Figures 6 et 7). 

Si Ton considere un objet 221 qui, dans une position 223, vient en 
interaction avec un autre objet 230, la direction privilegiee peut etre 
5 calculee uniquement sur des criteres geometriques (cas de la Figure 6) ou 
peut dependre de la configuration par laquelle les objets 223, 230 sont 
entres en collision en tenant compte d'un mouvement intermediaire 222 
d'au moins I'un des objets entre le pas de calcul precedent et le pas de 
calcul courant (cas de la Figure 7). 

10 Dans tous les cas le processus de detection des collisions permet 

d'extraire des couples d'elements d'objets en intersection, une longueur et 
une direction d'interference entre ces deux elements. 

Dans le cas privilegie d'une description de la surface des objets par 
des triangles, un element est soit un point, soit un segment, soit la face 

15 d'un triangle. La detection de collision peut prendre en compte trois cas 
canoniques d'intersection entre deux objets : intersection point/triangle, 
intersection segment/segment, intersection triangle/point. 

Le procede peut egalement fournir un ensemble d'elements d'objet 
en proximite qui pourrait potentiellement entrer en collision suite a 

20 Integration des forces de contact. Une distance et une direction separant 
ces elements sont alors calculees. 

Grace a la description de toutes les interferences et proximites 
entre les objets, I'ensemble des groupes de collisions de la scene peut etre 
construit. Chaque groupe va alors passer dans le module de contact. 

25 Apres la detection d'une collision, un traitement dans le module de 

contact permet de determiner la configuration du premier contact entre 
deux pas de temps de calcul. 

Si I'on prend en compte une description de la surface des objets par 
des triangles lorsqu'une zone d'interpenetration entre deux objets a ete 

30 definie a un instant T du pas d'echantillonnage de la simulation, il est 
extrait une liste de triangles composant la paire d'objets. Si I'on a un objet 
rigide et un objet deformable, les coordonnees du triangle representant 
Pobjet deformable sont traduites dans le cadre de reference de I'objet 
rigide a des instants distincts d'echantillonnage de simulation T et T-l. 

35 Pour toute paire possible triangle/triangle il est effectue une 

interpolation lineaire du deplacement de trois points Di, D 2 , D 3 du triangle 
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deformable entre les positions initiate et finale aux instants discrets 
d'echantillonnage T et T-l. On peut ensuite operer trois types de tests 
differents representes sur les Figures 3A a 3C sur la paire comprenant un 
triangle deformable 20 defini par des sommets Di, D 2 , D 3 et presentant 
5 diverses positions successives 21, 22, 23 entre les instants T-l et T et un 
triangle rigide 30 defini par des sommets Ri, R 2 , R3. 

Le test 1 (Figure 3A) correspond au cas ou un plan de collision est 
forme par le triangle rigide et va etablir une contrainte sur le point 
concerne du triangle deformable. 
10 Le test 2 (Figure 3B) correspond au cas ou un plan de collision est 

forme par un segment rigide et un segment deformable a I'instant de 
collision t, et va etablir une contrainte sur deux points de I'objet 
deformable. 

Le test 3 (Figure 3C) correspond au cas ou un plan de collision est 
15 forme par I'objet deformable a I'instant de collision et va etablir une 
contrainte sur trois points du triangle deformable. 

L'invention peut s'appliquer aussi bien au cas de collisions entre un 
objet rigide et un objet deformable qu'au cas de collisions entre deux 
objets deformables. 

20 Le module de contact va etre decrit ci-dessous de facon plus 

detaillee en reference au cas prefere d'une description des objets en 
triangle. 

Le module de contact est appele autant de fois qu'il y a de groupes 
de collisions dans la scene a I'instant du calcul. Le module de detection de 
25 collisions a stocke dans un espace memoire pour chaque collision : 

- la normale, 

- la paire d'objets et les elements concernes par la collision, 

- (eventuellement) les points d'application de la force de contact. 

Si I'algorithme de detection de collision ne donne pas les points 
30 d'application de la force de contact (cas illustre sur la Figure 6 d'une 
detection sans mouvement intermediate), il faut les reconstruire et dans 
tous les cas, il faut interpoler ces points d'applications par rapport au 
modele de deformation choisi. 

Dans le cas prefere d'une description des objets en triangle, on 
35 separe le probleme en deux cas : soit on a deux elements segments, soit 
on a un element nceud et un element face. 
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Dans le cas d'une intersection segment/segment entre un triangle 
40 defini par des sommets Pi, P 2/ P 3 et un triangle 50 defini par des 
sommets Qi, Q 2/ Q 3 (Figure 4), les deux points choisis 41, 51 sont situes a 
I'intersection de chacun des deux segments Qi, Q 2 , resp. Pi, P 2 avec le 
plan forme par la face de I'autre triangle en intersection. 

Le vecteur reliant les deux points trouves 41, 51 est appele 8. 

Dans le cas d'une intersection Point/Face entre un triangle 60 defini 
par des sommets Pi, P 2 , P 3 et un triangle 70 defini par des sommets Qi, 
Q 2 , Q3 (Figure 5), le premier point choisi est le point de I'intersection et le 
second est la projection de ce point sur la face de I'intersection. Le vecteur 
reliant les deux points trouves est appele 8. 

On peut utiliser I'algorithme d'intersection avec mouvement 
intermediate propose , par Xavier Provot (Collision and Self-collision 
handling in cloth model dedicated to design garment. Graphics Interface 
1997, 177-189) qui permet d'obtenir une configuration approchee entre 
les deux triangles au moment de la collision. 

Pour repartir les deplacements et les forces de ces points, dans le 
cas privilegie d'une interpolation lineaire pour la modelisation en elements 
finis (triangles - tetraedres), on utilise les coordonnees barycentriques. 

Pour pouvoir calculer correctement les deformations sur les objets, 
il faut garantir la non-interpenetration. L'interpolation des elements est 
done utilisee de fagon a n'avoir qu'une inconnue de force et de distance 
par contact. Pour simplifier le probleme, on se place dans le cas d'une 
interpolation lineaire pour les elements finis. 

On a ainsi dans le cas d'un contact segment/segment (Figure 4) la 
relation suivante : 



ou alpha et 1-alpha sont les coordonnees barycentriques sur le premier 
segment Qi, Q 2 et Beta, 1-Beta sur le deuxieme segment Pi, P 2/ 
ai, bi, q sont les coordonnees de la normale r\\ du triangle 40, 
Wi, W 2 sont les coordonnees du premier segment Qi, Q 2 , 
Vi, V 2 sont les coordonnees du deuxieme segment Pi, P 2 . 
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Dans le cas d'un contact point/plan (Figure 5), ai, bi , q sont les 
coordonnees de la normale rii du triangle 60, la distance d'interpenetration 
8 entre les triangles 60 et 70 s'ecrit : 
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ou : alpha, beta, gamma sont les coordonnees barycentriques sur le 
triangle 60 (addition = 1). 
10 Wi, W 2 , W 3 sont les coordonnees du premier triangle 60, et Vi sont 

les coordonnees du point de contact 71 constitue par un sommet Qi 
du second triangle 70. 

On utilise une interpolation identique pour la force de contact. 
Une fois trouve le point d'application des forces de contact, les 
15 caracteristiques mecaniques des objets sont transferees dans I'espace 
defini des contacts. Pour la suite, on suppose que Ton traite un groupe de 
m contacts avec n objets. 

Pour le transport des caracteristiques mecaniques globales dans le 
cas ou Ton considere la masse et I'inertie de I'objet de facon globale, en 
20 son centre de gravite, on peut utiliser une matrice jacobienne classique, 
definie dans les travaux de Ruspini (Diego Ruspini & Oussama "A 
Framework for Multi-Contact Multi-Body Dynamic Simulation and Haptic 
Display", Proceedings of the 2000 IEEE/RSJ International Conference on 
Intelligent Robots and Systems"). Grace a cette matrice, on a une relation 
25 instantanee entre les forces de contact f c dans la direction du contact, les 
accelerations 8" c (contraintes) et 8"iibres dans la m§me direction : 

<?> JeM-j/f^^ (3) 

30 ou J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane lineaire et angulaire dans I'espace des contacts, M est une 
matrice bloc diagonale correspondant a la masse et Tinertie des n objets 
du groupe de contacts et J C T est la matrice transposee de J c . 
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Pour le transport des caracteristiques locales, on peut utiliser 
I'interpolation definie pour les triangles. On a ainsi une relation entre les 
deplacements des points du maillage deformable et les deplacements dans 
I'espace des contacts et la meme interpolation peut etre utilisee entre les 
5 forces de contact et les forces sur les nceuds. 

En reprenant la relation Hneaire instantanee de caracterisation des 
deformations, on obtient : 

5 C = [S i=1 n N; A (U k .,) (N C ') T ] / c+ 8 libre (6) 
ou Nc represente la matrice de passage de I'espace des deplacements du 
10 maillage vers I'espace des deplacements aux contacts et A est une matrice 
comme defini plus haut. 

La modelisation du contact est choisie de facon a respecter au mieux 
les lois physiques. Par souci de simplification on peut neanmoins de 
preference considerer que les contacts ne changeront pas de direction 
15 durant la resolution du calcul meme si en pratique ce n'est pas strictement 
le cas. 

Le premier postulat du probleme de Signorini est qu'il n'y a pas 
d'interference entre les objets s'ils sont solides (les matieres ne se 
melangent pas). Ainsi, on souhaite qu'apres la resolution du probleme le 
20 deplacement au contact soit positif ou nul : 

S>0 (9) 

Le second postulat est que Ton est dans le cas d'un contact sans 
25 frottement, done la force de contact est dirigee selon la normale : 

/>0 (10) 

Le troisieme postulat est que la force de contact est non nulle (f c * 0) si et 
seulement si il y a reellement un contact {Si =0). Cela cree une relation 
30 complementaire entre les deux vecteurs : 

Sc-L/c (ID 



35 



Pour pouvoir resoudre le probleme de Signorini, on utilise le 
transport des caracteristiques mecaniques. Les effets des caracteristiques 
locales et globales sont additionnes en integrant les accelerations pour se 
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ramener uniquement a une relation entre emplacement et force dans 
I'espace des contacts. 

Avantageusement, on utilise une methode numerique qui a 
tendance a conserver I'energie comme la methode des trapezes (aussi 
5 nommee Tustin), pour integrer les accelerations et qui peut s'exprimer de 
la fagon suivante si Ton considere une grandeur X aux instants t (Xt) et 
t+1 (X t+ i) : 

X t+1 =X t + »/ 2 dt(X t ' + X t+1 ') 
io x t+1 *=x; + i/ 2 dt(x t " + x t+1 ") (12) 

En utilisant cette methode numerique et les equations (4) et (5), on 
obtient la relation suivante : 

15 S c = [i/4 dt 2 j c m- 1 j c T +s i=1 n n;c (u,,) (n;) t ] / c + \ bre (is) 

Le coefficient 1/4 peut etre different si Ton utilise une autre 
methode d'integration des accelerations. 

Si le modele mecanique choisi ne comporte pas de caracteristiques 
20 globales et si la masse et I'inertie sont integrees au niveau local dans le 
modele de deformation, on utilise uniquement I'equation (5) qui comporte 
deja implicitement une integration numerique en temps. 

S c = [E i=1 n N; A (U k ,) (N C ! ) T ] f c +8 libre (6) 

25 

Les postulats definis dans le probleme de Signorini et I'equation 
lineaire instantanee creant une relation lineaire entre forces de contact et 
deplacements dans I'espace des contacts permettent de formuler le 
contact sous forme de probleme de complementarite lineaire (LCP). Pour 

30 ce type de problemes, il existe de nombreux algorithmes de resolution 
(voir par exemple Murty, K.G., Linear Complementarity, Linear and 
Nonlinear Programming. Internet Edition 1997) qui sont capables de 
resoudre le probleme en un temps compatible avec les performances 
demandees par I'haptique. A titre d'exemple, on peut effectuer un calcul a 

35 une frequence de I'ordre de 500 Hz a 1000 Hz pour un nombre 
raisonnable de contacts (par exemple 30 a 40 contacts) en utilisant 
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ralgorithme du Pivot principal sur un ordinateur de type PC Pentium IV 
2GHz. 

La Figure 3 illustre 1'interaction entre un objet deformable 80, tel 
qu'une pince avec un autre objet 90. Dans ce cas, I'objet deformable 80 
5 est virtuellement attache a I'interface haptique (des lors qu'il est tenu par 
i'utilisateur) dans une zone d'un nceud O definissant un repere Ox 0 yoZo. 
Pour de tels points on peut appliquer les conditions de Dirichlet. 

A chaque pas de simulation, la configuration de mouvement libre 
donne un espace de contact. Une resolution LCP donne les forces de 
10 contact f, non nulles et un mouvement contraint en est deduit. Ces forces 
(illustrees par la normale U de coordonnees ai bj q sur la Figure 9) sont 
transportees au point O pour creer la force et le couple sur Tinterface 
haptique. 

Les Figures 10A a 10C montrent I'exemple d'un objet deformable 80 
15 constitue par un clip mis en place sur un tube 91. On voit differentes 
deformations du clip 80 dans des positions 81, 82, 83 differentes par 
rapport au tube 91. 
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REVENPICATIONS 

1. Precede de simulation interactive du contact entre au moins un 
5 premier objet deformable (201 ; 20 ; 40 ; 60 ; 80) et au moins un second 
objet (202 ; 30 ; 50 ; 70 ; 90 ; 91) avec un pas predetermine de temps 
d'echantillonnage d'un modele simule, 

caracterise en ce que : 

(a) on calcule au prealable les parametres decrivant les caracteristiques 
10 physiques de chacun des objets, tels que la geometrie et la 

mecanique des materiaux de chacun des objets, et on stocke ces 
parametres dans une memoire, 

(b) au debut de chaque pas de temps d'echantillonnage du modele 
simule, on procede au niveau de chaque objet a une analyse en 

15 temps reel du comportement propre de I'objet pour predire les 

positions, vitesses et accelerations de cet objet selon un mouvement 
libre qui ne tient pas compte d'eventuels contacts ulterieurs, 

(c) a chaque pas de temps d'echantillonnage du modele simule, on 
analyse en temps reel, au niveau d'une scene globale comportant les 

20 objets susceptibles de venir en contact, des paires d'objets qui sont 

detectes en intersection, et on etablit une liste de groupes de 
collisions qui contient une chaine d'objets en collision et une 
description des collisions, 

(d) a chaque pas de temps d'echantillonnage du modele simule, on 
25 rapatrie en temps reel, pour chaque groupe de collisions, des 

parametres representant les caracteristiques physiques des objets et 
la description des collisions, de maniere a determiner, pour chaque 
cas, la solution au probleme de Signorini qui gouverne le contact 
entre deux objets dans le cas d'un glissement relatif pur, 

30 (e) a la fin de chaque pas de temps d'echantillonnage du modele simule, 
on procede au niveau de chaque objet a une visualisation en temps 
reel du comportement propre de I'objet suite a la collision, et 
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(f) Pensemble des traitements en temps reel s'effectue avec un pas de 
temps de calcul plus court que le pas de temps d'echantillonnage du 
modele simule, de maniere a definir une simulation interactive ou 
I'utilisateur peut intervenir directement en cours de simulation. 

5 2. Procede selon la revendication 1, caracterise en ce que lors de 

I'etape a) de calcul prealable des parametres decrivant les caracteristiques 
physiques de chacun des objets, on utilise pour les parametres decrivant 
la mecanique des materiaux une description des deformations de type 
elements finis, avec le remplissage et I'inversion de matrices, la resolution 
10 de systemes d'equation et le stockage des donnees en memoire. 

3. Procede selon la revendication 1 ou 2, caracterise en ce que 
chaque objet est decrit dans une configuration au repos comme un 
ensemble de triangles reproduisant sa surface et un ensemble de 
tetraedres decrivant I'interieur de I'objet. 

15 4. Procede selon la revendication 3, caracterise en ce que chaque 

triangle est decrit par trois points, places dans un ordre qui permet de 
calculer des normales qui sont invariablement dirigees vers I'exterieur de 
I'objet. 

5. Procede selon la revendication 3 ou 4, caracterise en ce que les 
20 deformations des objets sont interpolees par la methode des elements 

finis en utilisant un maillage tetraedrique lineaire. 

6. Procede selon Tune quelconque des revendications 1 a 5, 
caracterise en ce que a chaque pas de temps de calcul, on integre lors de 
I'etape b) au niveau d'un objet, les forces explicites appliquees sur I'objet, 

25 qui sont deja connues au depart du pas de calcul, de maniere a definir le 
mouvement qu'elles creent sur I'objet, tandis que la valeur des forces de 
contact implicites, qui dependent elles-memes du mouvement des objets 
dans le pas de temps de calcul, est determinee lors de I'etape d) de 
recherche au niveau d'une scene globale, de la solution au probleme de 

30 Signorini. 

7. Procede selon I'une quelconque des revendications 1 a 6, 
caracterise en ce que lors de I'etape c) d'analyse au niveau d'une scene 
globale, on detecte geometriquement les intersections existantes entre les 
objets de la scene pour extraire des couples d'elements d'objets en 
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intersection, une longueur et une direction d'interpenetration entre les 
deux elements d'un couple d'elements d'objets. 

8. Procede selon la revendication 7, caracterise en ce que lors de 
I'etape c) d'analyse au niveau d'une scene globale, pour extraire des 
5 couples d'elements d'objets en intersection, une longueur et une direction 
d'interpenetration entre les deux elements d'un couple d'elements d'objet, 
on prend egalement en compte un mouvement intermediate des objets 
entre le pas de calcul precedent et le pas de calcul courant, pour calculer 
une direction d'interference privilegiee entre ces objets. 

10 9. Procede selon I'une quelconque des revendications 3 a 5 et selon 

la revendication 7, caracterise en ce que lors de I'etape d) de recherche de 
la solution au probleme de Signorini, on reconstruit les points extremes 
d'application de la force de contact entre deux objets soumis a une 
collision lorsque ces points extremes d'application n'ont pas ete 

15 determines a I'etape precedente. 

10. Procede selon la revendication 9, caracterise en ce que lors de 
I'etape d) de recherche de la solution au probleme de Signorini, dans le 
cas d'une intersection segment-segment de deux objets en triangle (40, 
50), les deux points (41, 51) choisis pour constituer les points extremes 

20 d'application de la force de contact entre les deux objets (40, 50) soumis a 
une collision sont situes a ('intersection de chacun des deux segments 
(P1P2/ Q1Q2) avec le plan forme par la face du triangle en intersection. 

11. Procede selon la revendication 9, caracterise en ce que lors de 
I'etape d) de recherche de la solution au probleme de Signorini, dans le 

25 cas d'une intersection point-face de deux objets en triangle (60, 70) , un 
premier point (71) choisi pour constituer un point extreme d'application de 
la force de contact entre les deux objets (60, 70) soumis a une collision 
est le point de I'intersection tandis que le second point extreme 
d'application de la force de contact entre les deux objets soumis a une 

30 collision est la projection (61) du premier point extreme (71) sur la face du 
triangle en intersection (60). 

12. Procede selon I'une quelconque des revendications 9 a 11, 
caracterise en ce qu'on utilise les coordonnees barycentriques pour 
repartir les deplacements et les forces des points d'application de la force 
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de contact entre les points extremes d'application de la force de contact 
en effectuant une interpolation lineaire pour une modelisation en elements 
finis. 

13. Precede selon les revendications 10 et 12, caracterise en ce que 
Ton calcule la distance 5 d'interpenetration entre les deux points extremes 
(41, 51) d'application de la force de contact dans le cas d'un contact 
segment-segment entre un premier segment (Qi Q 2 ) et un second 
segment (Pi P 2 ) d'un second triangle a partir de I'equation suivante : 



10 



8 = [a, b, cj 



[a I - a] 



w, 



- [p 1-P] 



(1) 



ou : a et 1-a sont les coordonnees barycentriques sur le premier 
segment (Qi Q 2 ), 

p et 1-p sont les coordonnees barycentriques sur le deuxieme 
15 segment (Pi P 2 ), 

ai bi q sont les coordonnees de la direction n\ d'interpenetration, 

Wi et W 2 sont les coordonnees du premier segment Qi Q 2 , 

Vi et V 2 sont les coordonnees du second segment Pi P 2 . 

14. Procede selon les revendications 11 et 12, caracterise en ce que 
20 Ton calcule la distance 8 d'interpenetration entre les deux points extremes 
(61, 71) d'application de la force de contact dans le cas d'un contact point- 
plan entre un point (71) d'un second triangle et un plan (Pi P 2 P 3 ) d'un 
premier triangle a partir de I'equation suivante : 







"w," 




5 -[a, b, cj 






- v, 






w 3 . 





ou : a, p et y sont les coordonnees barycentriques sur le premier triangle, 
ai bi q sont les coordonnees de la direction n, d'interpenetration, 
Wi, W 2 , W 3 sont les coordonnees du premier triangle (Pi P 2 P 3 ), 
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Vi represente les coordonnees du point de contact constitue par un 
sommet (Qi) du second triangle (Qi Q 2 Q 3 ). 

15. Procede selon I'une quelconque des revendications 1 a 14, 
caracterise en ce qu'apres avoir determine les points d'application des 

5 forces de contact entre deux objets en collision, lors de I'etape d) on 
procede au transfert des caracteristiques mecaniques des objets dans 
I'espace defini des contacts dans lequel on traite I'ensemble d'un groupe 
de m contacts avec n objets ournetn sont des entiers. 

16. Procede selon la revendication 15, caracterise en ce que lors de 
10 I'etape d) on considere la masse et I'inertie d'un objet de facon globale, en 

son centre de gravite et on etablit une relation instantanee entre les forces 
de contact f c dans la direction du contact, les accelerations S" c dues aux 
contraintes dans la meme direction et les accelerations libres 8" libre dans la 
meme direction connues lors de I'etape c) au niveau d'une scene globale, 
15 selon I'equation suivante : 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
20 instantane lineaire et angulaire dans I'espace des contacts, 

J I est la matrice transposee de J c , 

M est une matrice bloc diagonale correspondant a la masse et 
I'inertie de n objets du groupe de contacts. 

17. Procede selon I'une quelconque des revendications 9 a 14 et 
25 I'une quelconque des revendications 15 et 16, caracterise en ce que lors 

de I'etape d) pour le transport des caracteristiques mecaniques locales, on 
etablit une relation entre : 

• la difference de deplacement (U k ') des points du maillage 
deformable representant I'objet i a I'instant k, entre la 
30 deformation libre (U'uibre) et la deformation contrainte (U'k >c ) soit 

U k '= UV-UUre 
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• les positions relatives des objets, libres et contraintes, dans 
I'espace des contacts : 5nbre et 8 C , avec 

8 = Si=i n N c j U k ' + 5, ibre (4) 

5 

ou N c ' est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements aux contacts, et 

on etablit une relation entre les forces dans I'espace des contacts f c 
et les forces dans I'espace des forces de deformation F k :: 
10 F k = (N c i ) T /c (5) 

18. Procede selon la revendication 1, caracterise en ce que lors de 
I'etape d) on etablit une relation lineaire instantanee de caracterisation des 
deformations ou deplacement de contact 5 C a partir des forces de contact 
f c et des deplacements libres Snbre dus aux mouvements libres integrant 
15 uniquement les forces connues de facon explicite en debut de pas de 
temps de calcul, selon I'equation suivante : 

<? = [ V N e ' A (U,,) (N;) T ] f c+ % bre (6) 

ou : est une matrice de passage de I'espace des deplacements du 
20 maillage vers I'espace des deplacements de contacts, 

(N;) T est la matrice transposee de N^., 

A est une matrice permettant de definir la deformation de I'objet au 
niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et U k -! 
25 represente le vecteur de deplacement dans le repere local de I'objet 

au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 



U K = A (U k -0 F k + b (U k -!) (7) 

30 

ou F k est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, et 

b est un vecteur qui a une valeur dans I'espace des deplacements 
et qui depend du modele de deformation de I'objet. 
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19. Procede selon la revendication 17, caracterise en ce que lors de 
I'etape d) on etablit une relation instantanee de caracterisation des 
deformations ou deplacements de contact 8 C a partir des forces de contact 
f c et des deplacements libres Snbre dus aux mouvements libres integrant 
5 uniquement les forces connues de facon explicite en debut de pas de 
temps de calcul, selon I'equation suivante : 

S c = [e dr j c m 1 j/ + s i= ; n; a (u k ,) (N c y ] / c+ S libre (8) 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane lineaire et angulaire dans I'espace des contacts, 

y c est la matrice transposed de J c ' 

M est une matrice bloc diagonale correspondant a la masse et 
I'inertie des n objets du groupe de contacts, 

0 est une constante dependant de la methode d'integration en 
temps, 

N' c est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements de contacts, 

(N^) T est la matrice transposee de N;, 

A est une matrice permettant de definir la deformation de I'objet au 
niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et U k -i 
represente le vecteur de deplacement dans le repere local de I'objet 
au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 

Uk = A (U k -i) F k + b (U k -x) (7) 

ou F k est un vecteur representant les forces exterieures appliquees a 
30 I'objet exprimees dans le repere local, et 

b est un vecteur qui a une valeur dans I'espace des deplacements 
et qui depend du modele de deformation de I'objet. 
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20. Procede selon Tune quelconque des revendications 1 a 19, 
caracterise en ce qu'il comprend en outre une etape de couplage avec un 
module d'interface haptique pour produire un retour de sensation haptique 
sur un dispositif mecanique par lequel un operateur manipule les objets 

5 dans une scene virtuelle. 

21. Dispositif de simulation interactive du contact entre au moins 
un premier objet deformable (201 ; 20 ; 40 ; 60 ; 80) et au moins un 
second objet (202 : 30 ; 50 ; 70 ; 90 ; 91) avec un pas predetermine de 
temps d'echantillonnage d'un modele simule, caracterise en ce qu'il 

10 comprend : 

(a) un module (100) de calcul prealable des parametres decrivant les 
caracteristiques physiques de chacun des objets, tels que la 
geometrie et la mecanique des materiaux de chacun des objets, 

(b) une memoire (102) de stockage des parametres prealablement 
15 calcules dans le module (100) de calcul, 

(c) un module (101) de couplage avec une interface utilisateur (104) 
comprenant un dispositif mecanique tenu par un utilisateur lui 
permettant d'exercer virtuellement des forces sur lesdits objets dans 
une scene du modele simule, 

20 (d) un ecran de visualisation (107) pour visualiser lesdits objets 
representes sous forme de maillages, 

(e) une unite centrale de traitement (100) associee a des moyens 
d'entree (103), comprenant au moins 

el) un module d'analyse d'objet pour analyser en temps reel au 
25 niveau de chaque objet le comportement propre de I'objet pour 

predire les positions, vitesses et accelerations de cet objet selon 
un mouvement libre qui ne tient pas compte d'eventuels 
contacts ulterieurs, 

e2) un modele d'analyse d'une scene globale comportant les objets 
30 susceptibles de venir en contact, pour analyser en temps reel 

des paires d'objets qui sont detectes en interaction et etablir 
une liste de groupes de collisions qui contient une chaine 
d'objets en collision et une description des collisions, 
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e3) un module de rapatriement en temps reel, pour chaque groupe 
de collisions, des parametres representant les caracteristiques 
physiques des objets et la description des collisions pour 
determiner, pour chaque cas, la solution au probleme de 
5 Signorini qui gouverne le contact entre deux objets dans le cas 

d'un glissement relatif pur, 

e4) un module de traitement de chaque objet pour permettre en 
temps reel au niveau de chaque objet une visualisation en 
temps reel du comportement propre de I'objet suite a une 
10 collision, et 

e5) des moyens de determination d'un pas de calcul plus court que 
le pas de temps d'echantillonnage du modele simule de maniere 
a definir une simulation interactive. 

22. Dispositif selon la revendication 21, caracterise en ce qu'il 
15 comprend des moyens de production d'un retour de sensation haptique 

sur I'interface utilisateur (104). 

23. Dispositif selon la revendication 21 ou 22, caracterise en ce 
que le pas de calcul correspond a une frequence egale ou superieure a 
environ 500 Hz. 
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